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We present numerical results for electron spin relaxation rates for single and laterally coupled 
double GaAs quantum dots in a perpendicular magnetic field. As source of spin relaxation we 
consider hyperfine interaction with the nuclear spins in the GaAs substrate. Due to the differences 
in the energy scales of the nuclear and electronic Zeeman energies, the phonon bath system has to be 
taken into account for energy dissipation. The corresponding transition rates of second order show 
strong dependencies on correlations between the electrons and the electronic energy differences, and 
hence on the magnetic field. For a highly asymmetric double dot we have found a relatively low 
second order electron spin relaxation rate for a wide range of magnetic fields. 
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I. INTRODUCTION 

Strong efforts have been placed on testing and developing semiconductor quantum dots for quantum computing 
applications. For the implementation of a device that allows for operation on two spatially separated electron spins, 
i.e. a quantum gate, Burkard, Loss, and DiVincenzo have proposed two laterally coupled quantum dots containing 
one electron each.^ Such a quantum dot system can be designed with metallic gates on a hetero structure. With 
variation of the gate voltages a controlled, time-dependent inter-dot coupling can be achieved. This is an advantage 
compared to vertically coupled quantum dots. 

For (quantum) computation, the states representing the information should at least be stable during the calculation 
and the read-out afterwards. Hence it is very important and interesting to study the mechanisms of spin relaxation 
for these quantum dot systems. "^"^^ For magnetic fields of the order of IT and above, spin-orbit coupling is considered 
the most effective mechanism of electron spin relaxation in GaAs quantum dots,^'^'^ while for smaller magnetic fields 
the hyperfine interaction with the nuclear spins in the GaAs substrate is assumed to be dominant. -'^^ Johnson et 
al.^'^ have shown that up to decay times of 1ms, the triplet-singlet relaxation of two electron spins in a GaAs double 
quantum dot in a small perpendicular magnetic field is governed by hyperfine interaction with the nuclear spins. For 
a numerical study of relaxation of correlated electron spins due to spin orbit coupling see [20]. 

Erlingsson'^^ has calculated electron spin lifetimes for two non-interacting electrons in a single quantum dot due to 
a second order hyperfine/electron-phonon interaction process. We extend his approach to interacting electrons and 
to more complicated quantum dot structures, i.e. (asymmetric) laterally coupled double dots, allowing for a tayloring 
of electron spin relaxation rates. 

II. THEORY 

A. Transition rate formalism 

The electron spin relaxation rate due to interaction between the electrons in the quantum dot and the nuclear spin 
and phonon bath systems is calculated via Fermi's golden rule. At nonzero magnetic fields, the nuclear Zeeman energy 
is much smaller than the electronic Zeeman energy. Furthermore, for S"^ transitions also orbital energy differences 
need to be absorbed. Except for energetic crossings between the initial and final electronic states of the spin flip, 
the law of energy conservation forbids these "hyperfine-only" relaxations. The electron-phonon interaction in turn 
cannot directly attack the electronic spin degree of freedom. Yet, relaxation is possible in second order via a combined 
hyperfine and electron-phonon interaction process. Thermally averaging over the initial bath states, the transition 
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e( is the energy of the electronic state \£). Hu) is the energy of the phonon emitted during the relaxation process. 

The nuclear spin flip energy is neglected in the energy differences in Eq. (1), because it is much smaller than the 
corresponding electronic energies. P{n) and P{N) are the probabilities to initially find the nuclear system and the 
phonon bath system in the states In) and |iV), respectively. Eq. (1) implies that there initially are no correlations 
between the electronic subsystem and the bath systems. 

The initial state \i) and the final state |/) are the lowest energy states in the singlct/triplct sector. Eq. (1) diverges 
for the case of an energetic crossing between the initial state \i) and an intermediate state \£) with the same spin as the 
final state |/) (or the corresponding case for a crossing of final and intermediate state). The resulting vanishing energy 
denominators would need a thorough treatment. For the parameter ranges we studied, these types of crossings did 
not appear. Crossings between intermediate states do not cause pronomiced effects on the second order transition 
rate, because the contributions to the rate of the corresponding transition elements are continuously weighted with 
the energy differences to the initial and final state, respectively. 

The Markovian approximation involved in the usage of Fermi's golden rule generally is valid for the interaction with 
the phonon bath due to short phonon lifetimes^^ compared to typical electronic timescales in semiconductor quantum 
dots. It is known that the system of nuclear spins can be polarized by the hyperfine interaction with a 2DEG electron 
gas.^^'^^ In our case, however, the number of electrons is small. We can therefore assume that the nuclear system 
is not polarized for all times and the few dot electrons will not carry out enough collisions to polarize the nuclear 
system. 



respectively, r and Rk are the corresponding positions. For GaAs consisting of isotopes with 7 = 3/2 but diff'erent 

nuclear gyromagnetic ratios, we consider an average hyperfine constant A. From electron spin resonance data for 
conduction electrons in a GaAs layer of a GaAs/Ala;Gai_a;As heterostructure,^^ an approximate value for A of 
-2.5 • 10~^m^eV~^s~^ can be calculated. 

Projecting (2) on the electronic orbitals |n), we notice that the hyperfine interaction is modulated by the orbitals 
wavefunctions at the positions of the nuclei: 



The first sum, diagonal in the orbitals, only shifts the electronic eigenergies. We will take into account these diagonal 
terms by an effective magnetic field, which in addition to the external magnetic field shifts the electronic energies via 
a Zecman term.^°'^*'^^ The effective magnetic field for an electron with orbital quantum number n is calculated with 
the thermal meanvalue 



B. The hyperfine interaction 



For an s-band electron the hyperfine interaction is of Fermi contact type: 
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with the effective electronic Lande factor g*. Each Jfe has the same thermal expectation value for non-interacting 
nuclear spins: 
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where the magnetic field is pointing in z-directiou. 7 ^ ~ 58.5 • lO^s^^ T^-'^ is an average nuclear gyroniagnetic ratio. 
ks is the Boltzmann constant. Since the orbitals can be assumed to vary little on the nuclear distances, we can 
approximate the sum over the nuclei by a three-dimensional integral and have 

^eff«(/>T-I— -3' (6) 

with the lattice constant ao- We have dropped the index k since all nuclear spin expectation values are equal. 
Important for the validity of this semiclassical effective field arc a large number of nuclei and an unpolarized nuclear 
spin bath.^^ We notice, that Beff does not depend on the orbital quantum number of the electron. 

We treat the nondiagonal terms in Eq. (3) perturbatively, the electronic Zeeman energy scale being much smaller 
than the orbital energy differences. Introducing the creation and annihilation operators and a^^ for an electron 
in orbital n with spin quantum number a G {T, i}) we have 

VhF = E'( E 4a(^|5i^V') 

k nn'rra' \a£{x,y,z} 

x{n\Rk){Rk\n')a+^a^,„, . (7) 

With the prime on the sum over the electronic quantum numbers we denote that the sum runs over all quantum 
numbers except for which n = n' and a = a' . These diagonal terms are included in the effective magnetic field. The 
sum over the nuclei is again approximated by a three-dimensional integral. 

Evaluating the absolute square in Eq. (1), the hyperfine interaction operator (7) appears in products of matrix 
elements of the form 

E^(M)(o;M|VkF|6)(c|VkF|d;M) (8) 

(I a), 1 6), |c), and |d) are multi-electron states). Following Erlingsson,^* we redefine the nuclear spin bath state vector 
to only consist of the fluctuating pair /jt, 1^1 \ 

Im) ^ iMfe)® iMfc'), (9) 
assuming that the nuclear system always is unpolarized, since few electrons are not sufficient for a polarization. With 

^^a^{x,y,z) = la - {Ia)T (10) 

we write the fluctuations in the hyperfine interaction with the fcth nucleus due to nuclear spin fiuctuations as 

SV^^ = ^ E' E '54a(^l^i^V') 

{n\Rk){Rk\n')a+„a^,„,. (11) 
The sum over the initial nuclear spin bath states is replaced by a sum over all pairs of fiuctuating nuclear spins: 

E^(M)(a;M|VHF|6)(c|VHFM;M) 

^ EE^('"*=)('"'^l(«I^HF|&)(c|t>HFM)K)- (12) 
k Hk 
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P{^k) denotes the Boltzmann weight of a single (nuclear) spin in a static magnetic field (cf. Eq. (5)). Fluctuations 
among components of different nuclear spins are neglected, since the dipolar interaction energies are much smaller 
than the nuclear Zeeman energies. 

The correlations in Eq. (12) for different components of the same spin are given by (a, /3, 7 G {x,y, z}) 

{SiaSii3:^a)T = '^^al3iih{i^) j., (13) 
7 

for the same components by 

{dUL)., = {il)^ - {iSr- (14) 



C. The electron-phonon interaction 



The electron-phonon interaction enters the second order transition rate in form of the products 

V P(iV)(a; 7V|ypH|6)(c|VpHM; N) 5{t, - ef - hw) (15) 



(|a), |6), |c), and \d) are multi-electron states.) N enumerates the phonons in mode w. As we are studying electronic 
relaxation processes at low temperatures, we restrict the interaction to phonon emission. Hereby, the sum over the 
phonon bath states |A^) times the probability P{N) can be evaluated to a temperature dependend form factor 

F{T,u;) = [l-exp{-{kBT)-^hu}]-\ (16) 

The electron-phonon interaction Hamiltonian is 

VpH = j^^l^ Ku{q){a\exp{-iqr)\b)Q+^a+^a,,^ 

vaha 

+ h.c., (17) 

where is the creation operator of a phonon with wavcvector q and branch v (longitudinal (1) or transversal (t)). 
For GaAs, there are two important coupling mechanisms to acoustic phonons: piezoelectric (PE) and deformation 
potential (DP) coupling. While the former couples to longitudinal and transversal branches the latter at sufficiently 
low temperatures only couples to longitudinal branches. ^'^ The coupling coefficients therefore are 

«i(g) = cr^ + c°^ (18) 

K,{q) = C^- (19) 



The DP coupling coefficient is given by 
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Cr = {2nf'DJ^^ \q\'/\ (20) 



The deformation potential for bulk GaAs is £> = 6.7eV.^^ The mass density is pcaAs = 5.65gcm~^.^'^ The speed of 
sound for acoustical phonons for GaAs is ci = 5.11 • lO^cms"^ and Ct = 3.34- lO^cms"^.^^ The PE scattering depends 
inversely on the modulus of the phonon wavevector and therefore is dominant for small phonon energies:^^'^^ 



= {27rf'ehu^^^A.{0,cP) (21) 

with the PE tensor element hi4 = 1.2 • lO^Vcm"^.^^ The PE scattering depends strongly on the direction of the 

wavevector q relative to the coordinate system (expressed by the azimuthal angle and the polar angle 9, which is 
measured with respect to the layer normal). This anisotropy stems from the transformation of the PE tensor from 
the normal mode to the crystal coordinate system, resulting in the anisotropy functions^^ 

Ai = ZQa^0^^'^ (22) 

A, = 2 (a^/J^ + + aW) - UaO, 0). (23) 



FIG. 1: Sketch of a double dot confining potential in a-direction. The dotted lines define the intervals in which the potential 
is piecewisely defined (using splines). 



a = cos{<p) sin(^), f3 = sm{(f)) sin(^), and 7 = cos(^) are the direction cosines of the wavevector. 

The energy to be absorbed by the phonon bath determines the length of the wavevector for the longitudinal and 
transversal branches. The possible directions of the wavevector have to be integrated out. The energy E is absorbed 
by a phonon of wavevector q = Eh~^c~^. The corresponding phonon density of states is proportional to due to 
the surface of the sphere in q-space, that corresponds to the possible directions of the phonon emission. Eq. (15) 
further scales with a factor of q caused by the deformation potential coupling (20). For very small wavevectors 
q, the piezelectric coupling (21) is dominant. Therefore, the electron-phonon contribution to the transition rate is 
proportional to E for small energies and to E^ for larger energies (for very small energies of the order of fce^'j the 
form factor Eq. (16) has to be taken into account). 



Considering a magnetic field perpendicular to the layer plane, each electronic orbital is a product of a two- 
dimensional layer wavefunction and a one-dimensional wavefunction for the layer normal direction. For the calculation 
of the hyperfinc and electron-phonon interaction transition matrix elements, the electrons are assumed to be confined 
in the layer normal direction by an infinitely high potential well. This confinement is considered that much stronger 
than in the layer plane, that the groundstate of this direction can be assumed to be occupied only. 

The Hamiltion describing two-dimensional single-electron orbitals reads: 



with the effective electron mass m* (6.1 • 10~^^kg for GaAs^''), and the vector potential A (B = V x A). The two- 
dimensional layer orbitals are calculated in a discretized spatial basis, allowing for piecewisely defined (asymmetric) 
double dot potentials (see Fig. 1). In this way, the curvatures in the dot centers, their distance d, and the barrier 
height Eh can be defined independently. We define the dots to be parabolic near the dot centers. The confinement in 
y-direction is assumed to be parabolic. 

Outside the dot regions we consider the confining potential to be constant V{x,y) = Ec- We choose this constant 
so large that no unbound states are mixed by the Coulomb interaction to the low-energetic orbitals we consider. In 
this region we choose a sparser grid than in the dot region, because the wave functions decay slowly here. 

The Coulomb interaction between the electrons is calculated in the layer plane only, again requiring a strong 
confinement perpendicular to the layer. The corresponding matrix elements are calculated in momentum space, 
taking advantage of the inverse distance operator being diagonal in momentum space: 



Hereby, only two-dimensional integrals have to be evaluated instead of four-dinicnsional integrals in real space. 

Since the Coulomb interaction does not affect the electron spin, the resulting two-electron Hamiltonian can be 
diagonalized in blocks of fixed pairs {S'^,Sz)- The few-electron states are calculated by exact diagonalization in a 
Fock basis. 'Exact' means here, that all correlations due to the Coulomb interaction within the finite basis set are 



D. The electronic system 
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FIG. 2: Second order transition rate (a) from the spin down state to the spin up state of an electron occupying the lowest 
orbital of a GaAs quantum dot and logartlimic plots of electron-phonon and hyperfine transition matrix elements (b) (arbitrarily 
scaled). The fitted curve is a straight line with slope « 3. 



taken into account. The Zeeman energy corresponding to the two-electron spin component Sz simply adds to the 
orbital energies. 

We have considered single- and two-electron systems in single and double dots. The systems have a confinement 
of ?KJo ~ 2meV. The double dots have a distance of about 20nm between the dot centers, separated by an energy 
barrier of 20mcV (due to the energy barrier between the dots the double dot potential is deformed, such that the 
single-particle ground state is increased to hiOQ « 2.75meV). For the layer thickness of the dot system we have chosen 
a typical value of 2.8nm, corresponding to ten monolayers of GaAs.^^ All calculations have been performed for the 
temperature T = O.IK. 



III. NUMERICAL RESULTS 
A. Transitions in single-electron systems 



For single-electron systems, we consider Sz transitions of an electron in the lowest orbital. For B = Q the corre- 
sponding two states = ±1) are degenerate. The second order transition rate is zero here (see Fig. 2(a)), because 
the phonon density of states vanishes for q ^ Q. Furthermore, the operator cxp(iqf) approaches 1, such that different 
states are not coupled. At i? = a first order, "hyperfine-only" transition is possible. For arbitrarily small magnetic 
fields B > Q, the electronic Zeeman energy is always larger than the nuclear Zeeman energy, forbidding the first order 
relaxation due to the law of energy conservation. 

The second order transition rate increases with the magnetic field, since the electron-phonon and hyperfine transition 
elements increase. In addition, the energies of states with higher angular momentum decrease, leading to smaller 
energy dominators in Eq. (1). Therefore, these states contribute more to the transition rate with increasing magnetic 
field. 
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FIG. 3: Second order transition rates between the lowest energetic two-electron triplets with Sz = 1 and iS^ = in a single and 
a double GaAs quantum dot, respectively. 



Fig. 2(b) shows the electron-phonon transition element from the initial state to the state with the same total spin 
in the first excited orbital. The logarithmic plot reveals that the absolute square of the element is proportional to B^, 
i.e. to the cube of the electronic Zeeman energy. For small q, the electron-phonon matrix element is proportional to 
B due to dominating piezoelectric coupling. 

The corresponding hyperfine transition element (Fig. 2(b)) leads the system back from the intermediate state, 
which we have considered above for the electron-phonon interaction, back to the final state. The hyperfine transition 
elements grow mostly due to the fact that the magnetic field shrinks the electronic wavefunctions. The matrix elements 
are proportional to an integral over the product of four electronic wavefunctions (cf. Eq. (8)). For wavefunctions of a 
harmonic oscillator, such an integral is inversely proportional to the width of the wavefunctions. Since the width is 
reduced with increasing magnetic field, this explains the growth of the hyperfine transition matrix elements. 

The second order transition rate of a single electron in a double dot is similar to the single-dot case. We have found 
the rate to be two orders of magnitude larger than for the single dot. This is due to the energetic differences to the 
next intermediate states. In the double-dot case, this energy is smaller due to tunnel splitting, resulting in a larger 
transition rate. 



B. Transitions in two-electron systems 



Spin transitions in two-electron systems can occur within the triplet sector, i.e. between triplets with Sz G {0, 1}, or 
between triplet and singlet states. For a two (electron) spin system, Sz transitions in the triplet sector are decoupled 
from the orbital degrees of freedom. 



1. Sz transitions 



Fig. 3 shows the second order rates between the lowest energetic two-electron triplets with Sz = 1 and Sz = for a 
single dot and a double dot, respectively. Since the two triplets arc degenerate at i? = 0, the second order rate is zero 
there. Again, a first order transition is possible here, of course, which is not considered here. A transition within the 
triplet sector is a simple spinflip. As in the single-electron case it does not involve any orbital transitions. Therefore 
the rate has the same magnetic field dependence. 



2. transitions in single dots 

In contrast to the spin flip processes considered so far, the singlet-triplet transitions involve the orbital parts of the 
wave functions. 

Fig. 4 shows the second order S'2=l-triplet-singlet transition rate (a) and the energy spectrum (b) (cf. [39,40]) of 
two electrons in a GaAs quantum dot with "normal" Coulomb interaction (e = 13.1; in order to study the influence 
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FIG. 4: (Color online) Second order transition rate (a) between the lowest energetic two-electron singlet and triplet with Sz = 1 
in a GaAs quantum dot (normal Coulomb interaction: e = 13.1) and energy spectrum (b) (only the lowest energies are shown; 
for the rate calculations a few hundreds of intermediate states have to be included for convergence of the rates). 



of the interaction strength, we will rescale this value). At i? = the second order rate is finite, since the energies 
of the singlet and the triplet differ and thus, a nonzero phonon density of states enters the rate. In the range of 
the magnetic field we studied, the rate shows two minima and two maxima. This behavior of the rate is caused by 
the electron-phonon interaction, as the transition matrix elements have minima and maxima at the corresponding 
magnetic fields (see Fig. 5(a)). Though the triplet-singlet energy difference decreases for small, increasing magnetic 
fields, the electron-phonon transition matrix element grows with the phonon density of states for magnetic fields 
up to 0.5T. The reason for this is, that with the decreasing energy to be dissipated and the corresponding phonon 
wavevector of decreasing length, a geometric resonance with the dot size is approached. Hence, the corresponding 
single-particle transition matrix elements grow. For higher magnetic fields, approaching the triplet-singlet crossing, 
the vanishing phonon density of states dominates the transition rate. In the vicinity of the triplet-singlet crossings 
shown here, at B = 0.92T and B = 2.95T, it causes the second order transion rate to vanish (Meunier et al}'' 
have experimentally observed triplet-singlet relaxation caused by electron-phonon interaction with an enhanced spin 
lifetime near triplet-singlet crossings). Between the crossings, the rate follows the energy dependence of the gap 
between singlet and triplet. In the vicinity of the triplet-singlet crossing, resonances in the transition rate can be 
expected due to first order, "hyperfine-only" relaxation.''^ 

Due to the increasing hyperfine transition elements (see Fig. 5(b)) the second maximum of the second order transition 
rate is higher than the first maximum. The electron-phonon coupling thus mainly is responsible for the vanishing 
second order rate near the energy crossings, while the hyperfine interaction is responsible for the overall growth. The 
explanation for the growth of the hyperfine transition matrix elements is again the reduction of the width of the 
electronic wavefunctions with increasing magnetic field. 

At about 1.6T a crossing between the initial singlet state and the next singlet occurs. This leads to a pronounced 
effect on the transition rate: the rate "jumps" by about 75%. This is due to the drastic change of the character of 
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FIG. 5: Modulus square of the electron-phonon transition matrix element (a) from the initial electronic state to the state in 
the next orbital with equal spin (two electrons in a GaAs quantum dot; normal Coulomb interaction) and modulus square of 
the hyperfine transition matrix element (b) from the intermediate state to the final state. 



the initial state. At this magnetic field the cxcdtcd singlet state, from which the spin transition to the ground state 
is considered, changes. While the change of the electron-phonon transition elements at the corresponding magnetic 
field is not as dramatic (see Fig. 5(a)), the hyperfine transition elements show a "jump" there (sec Fig. 5(b)). This 
is amplified by a large electron-phonon rate which changes continuously at this transition. The crossing in the initial 
state does not lead to a resonance, caused by a vanishing energy denominator in the transition rate formula, as this 
intermediate state is not an intermediate state that allows for a second order transition to the triplet ground state. 

Wc have studied the behavior of the second order transition rate for systems differing in the strength of the Coulomb 
interaction (modeled by a variation of the dielectric constant e). The second order rate is minimal for the case of 
"normal" Coulomb interaction (e — 13.1). This is in clear distinction of what one would expect for first order rates. 
In the first order rates, only the hyperfine transition matrix elements enter. These are lowered due to correlations 
between the electrons caused by a stronger Coulomb interaction (e = 6.55). In the second order rate the correlation 
induced reduction of the hyperfine rate is compensated by the increasing electron-phonon matrix elements, which 
grow with the triplet-singlet energy difference. Therefore, an increase of the Coulomb energy need not necessarily 
lead to an increased lifetime of the electron spin, if the dominant relaxation mechanism is the combined hyperfine 
and electron-phonon scattering. 



3. transitions in symmetric double dots 

Fig. 6 shows the second order triplet-singlet transition rate (a) and the energy spectrum (b) of two electrons in a 
GaAs double quantum dot with "normal" Coulomb interaction. Since the lowest energetic singlet and triplets are 
separated from the next singlet and triplets, crossings of intermediate states with the initial state, if any, could only 
occur at very high magnetic fields, which we have not considered here. The "jump" we have found in the single-dot 



10 



b) 



> 




11 
10.5 
10 
9.5 

9 
8.5 

8 
7.5 

7 
6.5 



B[T] 









_„,--''*''''*'^ 1st singlet - 
— ""^ 1st triplet (82=1) -- ■ 




2nd singlet -- . 




2nd triplet (5^=1) 




3rd singlet 


1 


3rd triplet (S,=1) 



0.5 



1.5 



B[T] 



FIG. 6: (Color online) Second order transition rate (a) between the lowest energetic two-electron singlet and triplet with Sz 
in a GaAs double quantum dot (normal Coulomb interaction: e = 13.1) and energy spectrum (b). 



rates, will thus not appear in the following plots. 

As for the single dot, the triplet and singlet states are non-degenerate and thus the second order rate is finite 
for B = 0. Due to the Coulomb interaction the dots are essentially singly occupied. Thus the energy difference 
between singlet and triplet ground state is small, and thus, no geometric resonance can occur in the electron-phonon 
transition matrix elements. Therefore, the transition rate monotonously falls approaching the triplet-singlet crossing. 
Expectedly, the second order rate vanishes in the vicinity of the triplet-singlet crossing. Due to a nearly constant 
energy gap, the overall growth of the rate beyond the triplct-singlct crossing is due to the increasing hypcrfine 
transition matrix elements. 

Owing to the tunnel splitting, the energy differences between the initial and final states in the double-dot case 
are smaller than in the single-dot case. As this causes a smaller phonon density of states entering the second order 
transition rate, the rate for the double dot system is smaller than the transition rate for the single-dot system. For 
-6 = 0, e.g., the transition rate in the single dot-system is about one order of magnitude larger. Furthermore, the 
next singlet and triplet states are separated from the lowest ones, which means a lower contribution of these states 
to the second order rate. 

For weak Coulomb interaction (e — 39.3) the second order transition rate is similar to the case of normal Coulomb 
interaction. However, due to differences in the spectrum (see Fig. 7(a)) the rates at low magnetic fields are ^10 times 
larger than for normal Coulomb interaction. 

For strong Coulomb interaction (e — 6.55) the second order transition rate at small magnetic fields is about 10 times 
smaller than the rate in the case of normal Coulomb interaction, although the energy differences to the intermediate 
states are of the same magnitude (see Fig. 7(b)). There are two reasons for this. Firstly, the overall energetic distance 
between the initial and final state is lower for the case of strong Coulomb interaction. Secondly, correlations in 
the electronic wavefunctions due to the Coulomb interaction have the effect of a reduction of the transition matrix 
elements. 
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FIG. 7; (Color online) Energy spectrum the lowest energetic two-electron singlets and triplets with Sz 
quantum dot (weak Coulomb interaction: e = 39.3 (a) and strong Coulomb interaction: e — 6.55 (b)). 



1 in a GaAs double 



4- S transitions in asymmetric double dots 



Fig. 8 shows the second order two-electron triplet-singlet transition rate (a) and the energy spectrum (b) for a 
double dot with the right dot being energetically 2meV lower than the other dot. Due to this asymmetry, for B — 
(almost) no tunneling is possible between the left and right dot orbitals (we recall that the effective single-dot single- 
particle energy is about 2.75meV). This means a larger energy difference between the singlet and the triplet than in 
the case of a symmetric double dot, and thus, the second order rate is larger. 

For B — the two-electron ground state is a singlet that has a large electron density in the right dot, but also a 
non-negligible density in the energetically higher left dot. The lowest triplet shows a relatively symmetric occupation 
of both dots (see Fig. 9). This is plausible, since the first excited single-particle orbital of the asymmetric double dot 
for i? = mostly is given by the left dot ground state. Of course, the Coulomb energy of the triplet ((Vp) « 2.0meV) 
is much smaller than the Coulomb energy of the singlet ((Vb) ~ 4.5meV) due to the larger distance between the 
electrons. At a magnetic field of -Bcross — 0.84T the reduction of the Coulomb energy becomes that advantageous, 
that a triplet-singlet crossing occurs and the triplet becomes the ground state. With increasing magnetic field, the 
energy of the first excited single-particle state in the left and the right dot is lowered (angular momentum m = — 1). 
When the first excited state in the right dot approaches the ground state of the left dot, tunneling between these 
two states becomes possible. The tunnel-splitting leads to smaller energy differences between the initial and final 
electronic states. Therefore, the second order transition rate is reduced and remains relatively small for magnetic 
fields B> B 

cross ■ 

In Fig. 10 triplet-singlet transition rates for different degrees of asymmetry are compared. For small energy dif- 
ferences between the left and right dot, tunneling is also possible for small magnetics fields i? — > 0. This leads to 
small second order rates. For larger energy differences between the dots, only with increasing magnetic field tunneling 
becomes possible. 

Fitting the data of the rate (see Fig. 8) for small magnetic fields B < -Bcross, we find that the rate is proportional 
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FIG. 8: (Color online) Second order transition rate between the lowest energetic two-electron singlet and triplet with Sz = 1 
in a highly asymmetric GaAs double quantum dot (normal Coulomb interaction: e — 13.1) and energy spectrum. 
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FIG. 9: Electron densities of groundstate singlet and lowest energetic triplet with Sz = 1 for the highly asymmetric double 
quantum dot at B = 0. 



to Ag^,, where AsTr is the energy difference between the singlet and the triplet. Since the electron-phonon transition 
matrix elements entering the second order rate approximately depend on the cube of this energy difference, the result 
of the fit confirms the assumption, that the electron-phonon interaction is responsible for the decrease of the transition 
rate. 

Only in the vicinity of the crossing, the energy difference is small enough to allow for a first order, "hyperfine- 
only" relaxation. Hence, if the magnetic field of the crossing is avoided, the singlet and the triplet can be considered 
relatively stable within the corresponding range of the magnetic field. This tunable behavior is desirable for the 
manipulation of quantum gates. 
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FIG. 10: Second order transition rate between the lowest energetic two-electron singlet and triplet with Sz = 1 for different 
degrees of asymmetry. 



Our numerical results for electron spin relaxation rates in quantum dots reveal to be remarkably influenced by 
correlations in the electronic wavefunctions, caused by the Coulomb interaction between the electrons. The second 
order rates further strongly depend on the electronic energy spectra. 

While first order rates can be considered to be lowered due to correlations caused by a stronger Coulomb interaction, 
the influence on the second order relaxation rates, mediated by hyperfine plus electron-phonon interaction, is more 
complex due to dependence on energy differences between the electronic states. 

For energetic crossings between the initial and final electronic states, first order "hyperfinc-only" relaxation is 
possible. Therefore, resonances in the relaxation rate are to be expected for these crossings. In the very vicinity of 
such a crossing, electronic states with quite long lifetimes can be expected, if the energy difference between initial 
and final electronic state is too largo for a nuclear absorption (and the hyperfine plus electron-phonon interaction 
relaxation mechanism is dominant). 

Promising for the implementation of a quantum gate, reduced second order relaxation rates have been found for 
double quantum dots. These rates are much lower than for single dots because of the tunnel splitting. This small 
energy difference leads to a low phonon density of states entering the second order rate, and thus, to a low transition 
rate. For asymmetric double dots, this tunneling can be enabled with increasing magnetic fields, possibly leading to 
a relatively low second order rate for a wide range of magnetic fields. 
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